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We construct explicit solutions of a number of Stieltjes moment problems based on moments of the form Pj (») = 
(2m) ! and p 2 r) (n) = [(rn) !] 2 , r = 1 , 2, . . . , n = 0, 1 , 2, . . . , i. e. we find functions W$ (x) > satisfying J °° x" W$ (x)dx = 

(r) 

Pj 2\ n )- ^ is shown using criteria for uniqueness and non-uniqueness (Carleman, Krein, Berg, Pakes, Stoyanov) 

that for r > 1 both Pj r l(n) give rise to non-unique solutions. Examples of such solutions are constructed using the 
technique of the inverse Mellin transform supplemented by a Mellin convolution. We outline a general method of 

(r) (r) (r) 

generating non-unique solutions for moment problems generalizing pj 2(1). such as the product Pj («) • pi, («) and 



[(rn)\]P,p = 3,4,.... 
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1 Introduction 

This paper concerns solutions of the Stieltjes moment problem |jl|, ||], i.e. positive functions W(x) which 
satisfy the infinite set of equations 
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We previously met this problem when considering the properties of the so-called coherent states (CS) of 
Quantum Mechanics These CS, as well as their many generalizations, should satisfy the property of 
Resolution of Unity which is essentially equivalent to Eq. (p]) Standard CS lead to p(n) = n\ and 

W(x) — e~ x . Generalized CS lead to p(n)'s other than n\ [^and the solutions of Eq. ([]]), if any, must be 
studied in each individual case separately [^, [| 0, [l0|]. It quickly occurred to us that an efficient approach 
to this problem was to use the inverse Mellin transform method which allows one to establish many 
solutions of Eq. (jl|), either by analytic methods [ |Tl| ] or by extensive use of available tables [12, 13 1. As a 



byproduct of this method we have established that, for a large number of combinatorial sequences such as 
Bell and Catalan numbers, etc., the corresponding sequences p(n) are solutions of the moment problem 
Eq. (jl|) [ [Ti] , Likewise, sequences arising in theory of ordering of differential operators jlq | solve 
appropriate moment problems too. This admixture of quantum-mechanical, analytical and combinatorial 
features deserves a deeper study which we intend to pursue. 

Our paper is partly expository in character, and has the following structure: first we establish a link 
between the Mellin transform and the moment problem; next, in Section 3, we provide principal solutions 
to two moment problems, termed toy models. In Section 4 we discuss criteria for the uniqueness of 
solutions of the moment problem. Section 5 is devoted to the explicit construction of non-unique solutions 
of the toy-models. In Section 5 some generalizations of model sequences, together with their solutions, 
are reviewed. Section 6 is devoted to a discussion and conclusions. 

We dedicate this paper to Philippe Flajolet on the occasion of his 60th birthday. His pioneering ap- 
plications of Mellin transform asymptotics to the analysis of combinatorial structures / |7^ [Zijj , [7^ / have 
been a source of inspiration for us. 

2 The Mellin transform versus moment problem 

The Mellin transform of a function f(x) of the real variable x is defined for complex s by the following 
relation 

M \f(x);s] = J x s - 1 f(x)dx = f*(s) (2) 


( in this definition f*(s) is not the complex conjugate of f(s) ! ), and its inverse 5W — 1 is defined by 

M- l \f*(s);x} = ^. J f*(s)x- s ds. (3) 

C—ioo 

See Ref. [ ^0| for a discussion of the dependence of f(x) on the real constant c. Among the many relations 
satisfied by the Mellin transform we shall mainly use the following 

•M \x b f(ax h );s} = -aT^ f ( — ^ , a,h > 0. (4) 
h \ h J 

If M [f{x);s] = f*{s) and M [g(x);s] = g*{s), then 

00 00 

M-\f{s)g*{s)-, X ]=jf(^)g{t) d l=jg^)f{t) d l (5) 
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which is called the Mellin convolution property. Note that if in Eq. (g) both f(x) and g(x) are positive 

for x > then / f(j is a l so positive for x > 0. This means that the Mellin convolution preserves 
o 

positivity, an essential property when considering the moment problem. 
Eq.dl]), if rewritten for n = s — 1 as 



W(x) = M~ l \p(s- l);x], 



(6) 



is, if W(x) > 0, a solution of a Stieltjes moment problem. Thus according to Eq. (||) one can solve 
the Stieltjes moment problem by performing the inverse Mellin transform on the moment sequence and 
checking if the resulting function is positive. All the solutions in the sequel have been obtained using Eqs. 
@, (Q) and Note that W(x) obtained via Eq. (^) may not be the only solution of Eq.(|l]). We call 
W(x) > obtained via Eq. (^) from Eq. ([!]) the principal solution. 

3 Principal solutions of the moment problems 

Let us consider two sequences of integers given by 



p( r) (n) = (2m)!, n = 0,l,2,..., r=l,2,. 



(7) 
(8) 



p^(«) = [(r»)!] 2 , n = 0,l,2,..., r=l,2,.... 

In the following we shall obtain the solutions of the Stieltjes moment problem for the moment sequences 
given by Eqs. ([7]) and (^]), i.e. the functions w['\ (x) > satisfying 



j x n W^(x)dx^p i {l(n). 



(9) 



(r) 

From now on we shall refer to the model problems pj 2 {n) as toy models TM1 and TM2, respectively. 

a) TM1: we begin by obtaining w[ l \x). Note that (2n)! = T{2n + 1) = T{2{s - We now apply 
Eq.(0) with a = 1, b = — \ and h = i. We subsequently use M _1 [r(s);x] = e~ x which gives 



j x n w[ l \x)dx 










lyfx 



dx=(2n)\, n — 0,1,2, 



(10) 



In the same spirit we observe that (2m)! = T (2r (n + ^)) = T (2r (s — ^y))- Upon using Eq. (Q) but 



now with a = 1 , b = — ^4— ^ and h = J- one obtains 



J x n w[ r \x)dx = j x n 





1 


i " 











I2rx 2r 







dx={2rn)\, n = 0,1,2,. 



(ID 



We call (x) > the principal solution of the moment problem Eq. (fil] 
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b) TM2: we begin by deriving (x) — 94 _1 [r 2 (i);jc] and employ the Mellin convolution Eq. (||). 
By using the Sommerfeld representation of the modified Bessel function of second kind Kq(x) [ pT| ] one 
obtains 



J ' x n W^\x)dx = J x n 2K {2x?) dx={n\f, n = 0,1,2,.... (12) 
o o 

Subsequently note that [(r«)!] 2 = [r (r (n + \))] 2 = [r (r (s - ^j)] 2 and again apply Eq.(Q) with a = 1, 
ft = — and h = -. The result is 



x = / jc" 



j x"W^ r \x)d. 
o o 



^rM2x±-) 



dx= [(rn)!] 2 , n = 0,l,2,.... (13) 



Since Ab(f ) > for t > 0, W 2 (jc) > is the principal solution of the moment problem Eq. (pj[). 

4 Criteria of uniqueness and non-uniqueness of the Stieltjes mo- 
ment problem 

It was realized from the very beginning of the history of the moment problem that its solutions may not be 
unique; i.e. for a given moment sequence there may exist more than one solution. Stieltjes himself gave 
an example of a non-unique solution of a problem leading to what turned out to be the lognormal distri- 
bution [jffi]. This example, being about the only one available, was quoted repeatedly in the literature. As 
recently as twenty years ago new non-unique solutions of other types of problems have been constructed. 
Of special interest were Stieltjes moment problems arising in probability theory, related to investigation 
of probability distributions not determined by their moments. Consequently, the subject was further de- 
veloped and systematized, largely due to the comprehensive work of Berg p3[], Berg and Pedersen [^4j, 
Lin [p5||, Stoyanov [ p6[ p7| |, Pakes [|28|], Gut [^9| and others. From the practical point of view one needs 
criteria to decide whether the pursuit of non-unique solutions is reasonable. Such criteria are either based 
on the p(«)'s alone or on the solution W(x) alone, or on both p(«) and W(x), see below. From now on we 
assume that all the p(«)'s are finite and that W(x) is continuous. We now give , in a somewhat condensed 
form, a list of such criteria. 

CI Carleman uniqueness criterion (T. Carleman, 1922, [|l]j) This is based on the properties of the p(«)'s 
alone and is: 

IfS — Y<7=l [p( n )] = °°> men * ne solution is unique. 

This criterion does not imply that if S < °° then the solution is non-unique. In fact it is possible to 
construct models for which S <°° and solutions are still unique pq, pOj] . 



C2 Krein's non-uniqueness criterion (M. G. Krein, around 1950, [f3ll]) 

This is based entirely on the solution W(x) and does not involve the moments. 

V Jo ^i+^^ dx < oo, i.e., the so-called Krein integral exists, then the solution is non-unique. 
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C3 Converse Carleman criterion for non-uniqueness (A. Pakes, 2001, [28|, A. Gut, 2002, [|29[]) 
This is based on p(n)'s and W(x). 

If there exists x' >0 such that for x> x', 0<W(x) <°° and\y{y) = — In [lV(e- v )] is convexin (y 1 ,°°), 
wherey' = ln(^'), andif, inaddition, S = Y.™=i[p{ n )}~ lr ' <°°, then the solutionW (x) is non-unique. 
See |E3l E4L E6l E8l E9|l for various refinements of these criteria. 



5 Construction of non-unique solutions for TM1 and TM2 

(r) 

We first apply the criterion CI to sequences p j ,(n) (the logarithmic test of divergence of the series S is 

conclusive) and conclude that for r = 1 both solutions W^\x) — and W^'ix) = 2Ko(2x?) are unique. 

For r > 1 we find that £~ =1 [P j ^(n)] ^ is convergent. Application of criterion C2 gives convergence of 
the Krein integral, showing that the solutions of TM1 and TM2 are non-unique. These findings are 
confirmed by use of criterion C3: convexity of = — l n T \ (e x )] is proved as it is equivalent to 

e x l r > for x > 0, (TM1) and K { (2e x l r ) - K {2e x/r ) > for x > 0, (TM2). The quest for non-unique 
solutions for r > 1 is then well founded. 

We know of no general method to construct such solutions. However we have proposed a procedure 
based on the application of inverse Mellin transform which generates the required solutions [| 32 1 which 
we now briefly expose. 

The first step is to construct, within the framework of a given set of p(n)'s, a family of functions CDi. (x), 
parametrized by a constant k (to be defined below), such that all their moments vanish, i.e. 



Jx n &i l (x)dx — j x 6 1 &k(x)dx — 0, n = 0,1,2,..., s=l,2, 



(14) 



The Mellin transform of C0a(x), go£(s), vanishes for s — 1,2, Such functions are orthogonal to all 

polynomials and play an important role in the study of integral transforms [ f^ . For our purposes we 
choose a particular method of producing the functions CO^(jc): 



/ 



x n i£>k{x)dx = ^[\(n)-hk(n), (15) 



or equivalently 

J x s - l a k (x)dx = p[ r l(s-l)-h k (s-l), (16) 


where hk(s) is any holomorphic function vanishing for s = 1,2, Among an infinity of possible choices 

the simplest one is h k (s) = sin(nk(s + 1)) and it defines a discrete parameter fc = ±l,±2,±3, The 

function C0fc(jc) acquires new parameters now and is formally obtained by calculating the inverse Mellin 
transform 

ioo 
1 f 

®il,k( x ) = T~- / Pil(s- 1 )®n(nte)x~ s ds, fc = ±l,±2,±3,.... (17) 
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It turns our that for both Pj 2 (s) the integration in Eq.(|l7|) can be performed: 

a) TM1: for p\ r) {n) = (2m)!, (o[ r) k (x) is a special case of earlier evaluation ^ and reads: 



= -4 rre -^sin[fe7 t (V)+^ 1/2 ''tan(g)], (r>|Jt|) 
= W[ r) (jc) sin [kK ( 2^1 ) + x 1 / 2r tan ( § ) ] . 



(18) 



In Eq.(|T8|) we notice a pleasant factorization of (jc) 



b) TM2: for p^C") = [(rn)!] 2 the corresponding function d^l(x) has to be, in the first place, represented 

as 

(19) 



©2 i W = M ~ 1 [r(ra - 5+ l)r(rj - s + 1) sin (%ks);x\ 



which can be conceived as another case of Mellin convolution. A little thought gives the two partners to 
be convoluted as 



(20) 



(see Eq.( 18)) and 



II -» fa { $\x) = -±=re- x ' r sin [lac (==1) +x 1 / r tan (is) 



(21) 



(compare Eq.(|18|)). Thus, the integral form of (O^Jx) is 

«w=K /2) (f)«>g 2) (of 



(22) 



whose evaluation requires a number of changes of variables as well as the use of formula 2.5.37.2, p. 453 
of vol.1 Ref.[ 12 1, but is essentially elementary. The final result is 



-Re 



T (^ £ ^ko(2x 1 /2'-(i + /tan(f)) 



1/2 



(23) 



rx r 



where we note a "near" factorization of W 2 (jc). 

Armed with explicit forms for COj {(jc) and CO, Ax) we are in position now to write down the families of 
non-unique solutions. Their structure has the form: principal solution + const- G0&(jc). More precisely: 
TM1: 



w[ r \e,k,x)=w[ r \x) 



1 +esin laz 



2r-l 
2r 



x^ 2r tm I — 



(24) 



for real e, lei < 1. 
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TM2: 



forr>2|fc|. 



W^(y ) k,x)=wi r) (x) 



1+Y 



(25) 



As 



1+Y 



is an oscillating function of bounded variation, a constant y= j{k,r) can be always 



~ ir) 

found to assure the overall positivity of W 2 (y k,x). The above technique for obtaining non-unique solu- 
tions can be readily extended to moment sequences more general than pj r 2 (n). We shall simply mention 
two such extensions without entering into details. 

For pj (n) = [(rn)!] 3 we begin with the sequence (n!) 3 for which the solution is 



J x"MeijerG([ 



[0,0,0], [ }},x)dx=(n\f 



(26) 



where we use a convenient and self-explanatory notation for Meijer's G-function borrowed from that of 
computer algebra systems. The extension, via Eq. (Eh, leads to the principal solution 



fx" [_LMeijerG([[ ],[]], [[0,0,0], [ }},x^ 



dx = [(rn)!] 



(27) 



The integrand in Eq. ( 27 ) is a positive function for x > which cannot be represented by any other known 
special function. It possesses an infinite series representation in terms of polygamma functions, which we 
will not quote here. The Carleman sum S is convergent but the criterion C2 is not conclusive. Only the 
criterion C3 permits to ascertain the non-uniqueness. The corresponding function CO^Jt) is defined as 



co^(jc) = M~ l [r 2 (rs-s+ l)T(rs-s+ l)sin(rob);jt] 



(28) 



which can be calculated as Mellin convolution of I — > W 2 (x) and II — * tt> j r f ^ (x), see Eqs. (|l8|) and (p3[), 
respectively. 

As a final example consider the sequence P4 (n) = pj (n)pn (ra). The corresponding Mellin convolu- 
tion of two principal solutions 2 (•*•)' see ^1 s - ( H ) an d (H), yields directly the principal solution for 



fx" 




i-s/kx 



^ Er MeijerG([ 



dx 



= (2m)![(ra) 
n = 0,l,2..., 



M 2 



r=l,2,. 



(29) 



which is non-unique by C3 only, as C2 remains inconlusive. 
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6 Discussion and Conclusion 

We have demonstrated a methodology for obtaining unique and non-unique solutions of the Stieltjes 
moment problem using the Mellin convolution method. Although the initial moment sequences were 
simple and classical, one is rapidly forced to leave the realm of standard special functions, as the resulting 
solutions are special cases of Meijer G-functions. In most cases they resist the check for non-uniqueness 
via both the Carleman criterion CI and the Krein criterion C2, and criterion C3 appears to be the only tool 
to decide this question. We have generated parametrized families of non-unique solutions exemplified 

34[ 35|| and their properties 
are investigated. 



here by Eqs. dig) and (E3h. Such functions are now called Stieltjes classes | 
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